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1. Introduction 



Interest in the possibilities of large extra dimensions 0] and the solution to the 
hierarchy problem for non-factorisable geometries |0] , has sparked a renewed effort in 
Kaluza-Klein theories. In particular, much work has been done on one loop quantum 
effects in the Randall-Sundrum model [Q, ^, |, ^ ^, ^, with the hope of stabilizing 
the radion field. In [^] it was shown that massless bulk fermions lead to stability, 
but the parameters of the model have to be fine tuned to obtain the observed mass 
hierarchy. It was also argued that the radion mass was smaller than a TeV and 
ruled out by experiment. This was also confirmed for massive bulk fermions [0]. 
Consequently, classical or non-perturbative mechanisms have been saught to stablise 
the radion |1T0| . 

In fact, we shall argue below that the radion mass can be consistent with ex- 
periment if there is one additional fine tuned parameter in the model. This makes a 
total of three parameters to be fine tuned, as opposed to just two in this context (the 
cosmological constant and the TeV Higgs mass) in the standard model. However, it 
has long been a mystery why the effects of quantum gravity do not introduce Planck 
scale corrections into the standard model. 



In the earlier work on Kaluza Klein theories 11, 12, uM it was found that the 



kinetic terms in the effective action played an important role. Although they have 
no effect on the size of the internal dimension, they can lead to instability in what 
would otherwise be a stable compactification. 

In this paper we calculate the finite one loop correction to the radion kinetic 
term from massless fermions and find that it vanishes to second order in derivatives. 
The calculation adapts a diagramatic method for obtaining derivative expansions of 
a heat kernel 0, which is closely related to the 'covariant perturbation theory' 



first devised by Vilkovisky |T^, |I8|, [T^ . In this case the fine tunings are protected 
at one loop order. 



2. The radion action 

We will give a simplified picture of how one obtains the radion field from the five 
dimensional action (along the lines of [^). The fifth dimension will be taken to be 



an orbifold S^/Z2, then the full set of classical field equations can be obtained from 
an action 



^ = -^ / (^-2A)--^ / K- C^, (2.1) 
SttGs Jm 47rG5 JdM JdM 

where the boundary dAi consists of a hidden brane S/j and a visible brane S^, with 
extrinsic curvature scalars K. The Lagrangian density Cm represents brane matter 



fields, and the action correctly reproduces the brane boundary conditions |21|. We 



shall be mostly concerned with the vacuum energies Vh and V^. In the original 
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Randall-Sundrum model, Vh = —V^ = Vq, where 
and A = — 6k^. 

We will identify the radion with the relative motion of the branes on a fixed 
Randal- Sundrum background 

ds'^ = aivYg^^dx^dx"" + rfy^ (2.3) 

with a(?/) = exp(— Ky) and Ricci scalar R = —20k^. This assumes that the brane 
motion and the five dimensional gravity waves decouple at leading order. 

The moving branes are located at y = yh{x) and y = yv{x), with extrinsic 
curvatures given by Kh = —K^, 

Kv\g\^^'^ ~ —4:Ka'l — ay d'^a^ — Ka^^^da^f (2.4) 

where is the d'Alembertian operator. Substitution into the original action gives 
a reduced Lagrangian C = T — V , where 

T = \Cy>^{dayf - \Cht^\dahf (2.5) 

V = ByK^at - Bhti^ai (2.6) 

and 

ByK^ = Vo + V. (2.7) 

BnK^ = Vo - V, (2.8) 

CyK^ = 2Vo + V. (2.9) 

ChK^ = 2Vo - Vh (2.10) 

Provided that the motions are small, the reduced action will generate a consistent 
set of field equations. The equations are incomplete because the metric variations 
have been restricted, but for spatially homogenous fields the field equations can be 



completed with a single constraint T + V = |20 



The relative motion of the two branes can be isolated by introducing a change 
of variables, 

= yv-Vh (2.11) 

a = C~'(Chal-CXy^' (2.12) 

where C = {ChCvY/"^. In the new variables, 

T = \ZTK\daf - \CK-\daf (2.13) 

V = Zv (Bh - B,e-*^'^) (2.14) 
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where 

Zt = 2C^a\Che'"' - C^e-'"')-^ (2.15) 
Zv = C^a^e^^'^iChe^'' - C^e-^'')-^ (2.16) 

The negative kinetic term is associated with the fact that a represents a gravitational 
degree of freedom of the double brane system, and corresponds to a Friedmann 
equation of the usual form. The Friedmann equation in a„, which still determines the 
expansion rate of the visible brane, has non-standard signs leading to doubts about 



the consistency of the Randall- Sundrum model and standard cosmology |£2[. The 
crucial idea here is that when the two branes are tied together by Casimir forces, the 
dynamical equations are simplest when expressed in terms of the collective coordinate 
a and predict the usual cosmological evolution. 

The classical theory has the shortcoming that the potential does not have a 
minimum. This has lead to the consideration of one loop effects [Q. Calculations of 
the vacuum energy for gj massless fermion fields leads to a potential 

V = Zv (Bh - 5^e-^"" + Ae-^^^il - e"'^")-^) (2.17) 



A=^-h^. (2.18) 



where A is given in terms of the Riemann zeta function, 

30%/ 
1287r2 

An equilibrium configuration requires both V = Q and V = 0. This implies that 
> A and gives a relationship between B^ and B^ 

The ratio of the mass scale and the Planck scale on the brane is set by A = e"*"^. 
In order to obtain a mass hierarchy close to the measured value we need the value of 
B^ to lie very close to the value of A. The mass of the radion (to a relative accuracy 
of A) is then 

ml = AC-^K-^X" (2.19) 

As Garriga et al. have pointed out, a radion with this mass should have been 
observed in particle experiments before now. However, at the expense of one further 
fine tuned parameter, the radion can be invisible with a large mass for small values 
of Cj;, which corresponds to selected ranges of Vq. A consistent radion model requires 
three fine tuned parameters, 

^ -2k^A (2.20) 
Vh ^ K.'^A (2.21) 
Vo ^ (2.22) 

We will not attempt to explain why the Lagrangian should have such fine tuned 
parameters. However, the importance of the kinetic terms in fixing the radion mass 
has lead us to examine the one loop corrections to these terms in order to test the 
robustness of the fine tunings at the one loop level. 
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3. The Basic Method 



We will obtain a perturbative expansion of the one loop effective action for massless 
fermions on the moving brane background. The coordinate system is chosen so that 
the coordinate in the fifth dimension is constant on the boundaries, and then an ex- 
pansion in derivatives of the metric is performed. An analysis of the different fermion 
boundary conditions can be found in [0]. We shall consider the case of fermion com- 
ponents which satisfy Neumann boundary conditions, but the other possibilities can 
be treated in an equivalent way. 

The conformal invariance of the massless fermions allows us to simplify the prob- 
lem by using the conformally related fiat background metric, 

ds"^ = dx^ + dr^ . (3.1) 

with Th{x) < r < We can replace r with the coordinate 6 which is constant 

on the boundaries, 

e = 7r{T- Thix))//3ix) (3.2) 
where (3{x) = t^{x) — Th{x). The metric becomes 

ds^ = dx^ + p^n-^de + A^dx^'f, (3.3) 



where 



for = ^ - f and 



= 9,a^ + f 6^ (3.4) 



= /3-\d,/3) (3.5) 
bf, = p-'d^{Th + r,) (3.6) 

The unperturbed configuration is therefore A = and f3 constant. 

The one loop effective action W can be related to the Laplacian A of the metric 
( |3.3|) . We shall obtain a derivative expansion of W using heat kernel techniques. The 
heat kernel K{x, 9, x', 9', t) satisfies the equation 

dK 

^~AK = 6{x- x')6{9 - 9')6{t) (3.7) 

The explicit form of the Laplacian is 

A = -np-\d, - deA,)P7r'\d^ - A^de) - rr^p-'d'e = Aq + A,. (3.8) 

where Aj denotes terms depending on derivatives of (3, Ai with one derivative, A2 
with two derivatives etc. The first terms are 

Ao = -dl - Tc'p-'dl (3.9) 
Ai = -e^7r2(9^/?-2)92 + A^ded'^ + A'^d.de, (3.10) 
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where ^ = x' — x. 

It is possible to show using perturbation theory that the heat kernel has a 
representation in the form of a time ordered exponential. In bra and ket notation, 



K{x,e,x',9',t) = (x'9 



-Aot rp 



exp 



Jo 



x9 



(3.11) 



The calculation of the heat kernel can be performed most easily in momentum space, 
therefore we introduce momentum basis states \kn) and vertex operators 



Vmnhk' = (A;'m|e^°*'A,e-^°*>n). 
For Neumann boundary conditions, 

K{x, 9, X, 9, t) = [ V cos m9 cos n9 e^^-*-'^'* (t exp - / V{t')dt' 

J V'T^) m,n V J 



(3.12) 



(3.13) 



The integrated trace of the heat kernel 



d X 



d k ^ — ^ 



T exp- [ V{t')dt') . (3.14) 



As in reference [0], we arrange the time ordered exponential in such a way that we 
can use Wick's theorem. Defining 

d'^k 



(...) 



Zit)J (27r)\^ 



with Z{t) = (Ant)-^ Ene""^"*, gives 



K{t) = J d^xZit) (Texp (^_e-'='*/Ve'='*/2 



(3.15) 



(3.16) 



In order to construct the vertex operator we take the previous expresions for A/ and 
make the replacements 



9r 



9, dl 



where 6'(t') is defined in appendix B, b'^ = —2ik^ and 



-n 



._d_ 

'dk„ 



2ikn' + ik^t. 



(3.17) 



(3.18) 



The expansion of the time ordered products leads to contractions between the above 
operators which can be interpreted as propagators (see appendix A). 

To begin with, consider Th + = 0, i.e. = 9ca^. If is substituted into 
( p.l2| ) one finds the vertex operators 



(3.19) 
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-H- 




Figure 3: The non-zero diagrams: (a) gives a total divergence; (b) and (c) cancel out. 



where u;^ = vr^n^//?^. 

The last term in the expression for V2 can be removed by the vertex ehmina- 
tion trick described in [Q. This allows terms to be eliminated by modifying the 
contraction 

B (ti, tj) = 2S{U - tj) - 2t-\ (3.21) 

The remaining vertex V2, and all higher order vertices, can be obtained by differen- 
tiating Vi. 

veiJje^eDperators can — 
reBented jdiagramatically by wticesj 
wick-ex<CTiiaWegs representingHh£--(9^7^^^ 
operators (see figure |l| and figure |^). 
Each diagram contributes to the heat 




3^e: 



^igure'2: Diagrams corresponding tc/the/vertex 
operator V2. Each additional line represents an 
Wimm^^:^ fflapte corresponding to the two or equivalently to the 

. ^, ^ ^ ^, zeta function via a Mellin transform, 

terms m the vertex operator Vi. ihe Ime repre- 

sents the 6 operator and^l^^^li^e jw-^A^ijo^-iz^t^^^dt. (3.22) 
represents the 6' operator. J r(s) J 

The one loop effective action is given hj W = |C'(0)- 



4. Second order calculation 

We now perform the calculation of the contribution to the kinetic terms to second 
order. There are only three diagrams that contribute to the effective action. Other 
diagrams are zero due to the proper time integrals (see appendix A) or theta integrals 
(see appendix B). 

The first diagram is shown in figure ^ with one second order vertex and one 
contraction D{ti,ti). This diagram gives a contribution to the heat kernel of (see 
equation ( p.l6D ), 

Z{t){{3a)) = -lj^2(9'a,-2a')h'Y.e'^M^ (4-1) 

where the first factor is a symmetry factor and the first integral in (A7) has been 
used. The contribution of the diagram to the zeta function (denoted by (a) is given 
by equation (|3.22| ), 

Ca{s) = / d'x^id^a, - 2a')Y^ul--'\ (4.2) 
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The sum can be expressed in terms of the Riemann zeta function and the derivative 
of the zeta function at s = is then 

C(0) = -Y^yC;?(-2) / d'xp-'id'^a, - 2a'). (4.3) 
Using = (3-'d^(3 and ('^{-2) = -Cii(3)/(47r2), we get 

C(0) = -I^^Ci?,(3)vr^/ dHd'r'), (4.4) 

This term is a total divergence and makes no contribution unless the 3-branes have 
a boundary. 

The second term corresponds to the diagram in figure ^b. Here identities in 
appendix A and B have to be used. The diagram gives 

Z{t)m) = j^,d^a,h'Y.e--'u;l (4.5) 
The contribution to the effective action is now 

cm = Y^f^C;^(-2) / d^xf3-'{d^a,) 

1 Air'^ r 

-a-2) d'xp-\dpf, (4.6) 
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where we have used equation (|3.4| ) and integrated by parts. 

The third term corresponds to the diagram in figure For this diagram 

Z{t){m = ~lj^,a'^t'Y.e-M- (4.7) 
The contribution to the effective action is 

C(0) = -^^Q-2)/rf^x/3-V 

1 A 2 

-C'^{-2) [ d'x(3-\d(3f. (4.8) 
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Thus we see that equations ( [4.6|) and ( [4.8|) cancel. A similar calculation shows that 
the terms involving also cancel. 



5. Conclusion 

The method which we have used to evaluate the kinetic terms in the radion action 
for a brane world model can be extended quite easily. We have explicitly shown 
that there are no one loop second order corrections from massless fermions. It also 
follows by a simple extension of the arguments that there are no one loop second 
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order corrections from any conformally invariant bulk fields on a conformally flat 
background or from massive fields on a flat background. The results are therefore 
not restricted to the Randall-Sundrum model and could be used for other brane 
world scenarios. 

We might ask whether there are any derivative terms in the one loop corrections 
to the radion action. In fact, they do have to exist because of the renormalisation 
scale dependence. The renormalisation scale enters along with the 05 heat kernel 
coefficient |2^, which is of order K'^ in the extrinsic curvature and therefore eighth 
order in derivatives. We do not know yet whether this is the leading term in the 
derivative expansion. 

Although the calculation of one loop effects may have a bearing on the hierar- 
chy problem, the quantum corrections can also be considered in contexts which are 
independent of the hierarchy problem. For example, the cosmological evolution of 
these models appears to offer advantages over other stabalisation mechanisms in the 



way that the expansion couples to the energy density on the brane |Q. It would be 
interesting to persue this point in more detail. 

A. Propagators 

There are three different propagators corresponding to the distinct combinations of 
the operators, 

{T5^{UW{t2)) =5^"'D{UM) (A.l) 

{T5'^{t,)5^{t,))=5,^^{t,M) (A.2) 

{T5'^{hmh))=5,,D{hM). (A.3) 



These can be evaluated using the creation and annihilation operators defined in [|14 
leading to 

ta) = 2min(ti, ts) - 2t-\t2 (A.4) 
t2) = 2e{t2 - ti) - 2t-H2 (A.5) 

where 6(t) = 1 for t > and zero otherwise. The third contraction evaluates to 
—2t^^, but we replace ( [A.3| ) with 

D ih,t2) = 25{ti-t2) -2t-^ (A.6) 

in order to reduce the number of vertices. The relevant time integrals of the propa- 
gators are: 

r D{h,h)dh = t D{h, t2)dtidt2 = 

Jo 6 Jo o 
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I5(ti,t2)dtidt2 = 

Jo 

[ 0ih,ti) + ^ih,ti))dh = O 
Jo 

I D {ti,t2)dtidt2 = 
Jo 

t ^ 1 

Max{ti,t2) D iti,t2)dtidt2 = --t^ 

i 



B. Theta integrals 



f\il5(ti,t2)dtidt2 = 
Jo 

f h0{ti, h) + h))dh = -\t' 
Jo o 

{ti+t2) D {ti,t2)dtidt2 = 



t ^ I 

tit2 D iti,t2)dtidt2 = -t^ (A.7) 
6 



The 6 direction is special because the modes in this coordinate are discrete, {0\n) = 
\/2cos{n9) for Neumann boundary conditions and {9\n) = \/2sm{n9) for Dirich- 
let boundary conditions. We need a way to replace a 9 sandwiched between e^°* 
terms inside the \kn) matrix elements. For this we use the Cambell-Baker-Hausdorff 
formula, 

e^«*^ee-^°* =^(t). (B.l) 

where 9(t) is an operator 

9{t) = 9c- 271^ (3-Hde (B.2) 
An important contraction identity we require for figure 3c is 

{n\T9{ti)de9{t2)de\n) = ~^f3W^ + 4cj>ax(ti, ts) ± ^'(ti + ^2) + 4cu^tit2. (B.3) 

The results for other 9 dependent terms in the matrix elements are given by 



{n\n) = 
{n\9cde\n) = ± 

{n\{9,d0)'\n) = 



{n\de\n) = 



12 



{n\9l\n) 
{n\dg\n) 



TT 
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(B.4) 



The upper signs correspond to Neumann boundary conditions and the lower to 
Dirichlet boundary conditions. 
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